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ABSTRACT 


The  Naval  Academy  Heat  Balanced  Engine  (NAHBE)  has  been  studied 
extensively  for  the  last  five  years.  Efforts  to  understand  and  predict  the 
process  occurring  in  the  engine  were  complicated  and  difficult  to  verify 
experimentally.  Optimization  of  the  output  of  the  engine,  though  possible, 
took  much  too  long  a  time  through  parametric  variation.  An  answer  to  the 
problem  of  predicting  the  time  dependent  reaction  occurring  in  the  engine 
was  sought.  Also,  some  easy  design  tool  for  optimizing  the  reaction  was 
required.  One  possible  solution  to  these  problems  was  to  model  the  NAHBE 
piston  cap  through  the  Helmholtz  Resonator  Theory.  This  theory  generalizes 
the  time  dependent  reaction  of  a  volume  of  fluid  in  a  chamber  by  the  acoustic 
analog  to  the  mechanical  system  forced  oscillator.  The  theory  enables  one  to 
predict  or  design  the  acoustic  impedance  of  a  cavity  and  its  natural  frequency 
by  way  of  the  cavity  and  inlet  geometry.  This  study  shows  how  the  solution 
to  the  Helmholtz  Theory  may  be  obtained  and  the  fluid  mass  displacement  of 
NAHBE  cap  modeled  by  the  Helmholtz  Theory  predicted.  The  solution  of  the 
NAHBE  cap  reaction  by  the  Helmholtz  Theory,  however,  is  only  a  hypothesis  until 
it  is  substantiated.  Conclusive  evidence  is  shown  in  this  study  through  a 
high-speed  Schlieren  film  analysis  of  the  actual  combustion  process  in  a  NAHBE, 
that  the  theory  holds  true  for  and  can  predict  reality.  Thus,  it  is  proven 
that  the  Helmholtz  Theory,  which  can  be  used  as  a  design  tool,  can  model  the 
real  systems  found  in  heat  balanced  engines  today. 
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NOMENCLATURE 


Combustion  cylinder  cross-sectional  area 
Opening  radius.  Speed  of  Sound 
Damping  constant 

Acoustic  capacitance,  capacitance 
Speed  of  sound 
Gravitational  constant 
Combustion  cylinder  length 
Current 

Runge-Kutta  constant 
Ratio  of  specific  heats 
Cap  thickness,  neck  length 
Effective  length 
Runge-Kutta  constant 
Mass  volume,  mass  element 
Pressure 
Quality  number 

Acoustic  resistance,  resistance 
Opening  cross-sectional  area 
Temperature 
Time 

Volume  velocity 
Volume,  voltage 

Volume  displacement,  displacement 

Variable 

Variable 


Acoustic  Impedance 
1 


z 


displacement 
Compression  ratio 
G  Non-dimensional  distance 

p0  Density 

t  Non-dimensional  time 

to  Frequency 
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INTRODUCTION 


Growing  concern  over  rising  petroleum  fuel  prices  has  increased  the 
interest  in  economical  consumption.  To  solve  the  long  term  problem  of 
limited  resources  it  will  be  necessary  to  do  away  with  common  internal  com¬ 
bustion  engines  and  all  other  power  sources  requiring  fossil  fuels.  In  the 
interim,  however,  it  will  be  necessary  to  stretch  out  what  reserves  we  do  have 
so  that  technology  and  science  will  have  time  to  find  a  solution  to  the 
problem. 

Little  had  been  done  before  the  1970's  to  increase  fuel  economy  except 
in  the  field  of  aviation  where  extra  fuel  meant  less  payload.  The  basic 
design  and  concept  was  still  never  improved,  however.  Two  thermodynamic  models 
for  internal  combustion  cycles  that  have  been  known  for  some  time  are  the  air 
standard  Otto  and  Diesel  cycles.  The  Otto  cycle  is  a  process  which  uses  a 
heat  addition  at  constant  volume  to  result  in  work.  The  Diesel  cycle,  on  the 
other  hand,  is  a  constant  pressure  heat  addition  process.  A  combination  of 
these  two  cycles  can  be  done  in  equilibrium  or  quasi-equilibrium.  The  former 
is  known  as  the  classic  Dual  cycle.  The  quasi -equilibrium  cycle  was  found  to 
be  desirable  because  of  theoretical  increases  in  fuel  efficiency  and  decreases 
in  peak  pressures  (Reference  1).  The  quasi-equi librium  cycle,  however,  must 
meet  certain  limitations,  making  it  seem  impractical  at  first  glance. 

The  thermodynamic  model  for  the  quasi-equilibrium  process,  known  as  the 
Heat  Balanced  cycle,  is  summarized  in  figure  1  and  table  I  (Reference  2).  Two 
points  should  be  noted.  Firstly,  the  total  mass  is  separated  into  two  regions. 
One  region  has  heat  added  in  a  constant  volume  phase  (Otto  cycle)  and  the  other 
has  heat  addition  during  a  constant  pressure  phase  (Diesel  cycle).  Secondly, 
the  two  regions  must  be  kept  separate  and  nonmixing  if  quasi -equilibrium  is  to 
exist. 
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The  Naval  Academy  Heat  Balanced  Engine  (NAHBE)  obtains  the  desired 
results  by  placing  a  cap  on  a  normal  piston  so  as  to  create  a  chamber  of  gas 
partially  separated  from  the  combustion  chamber  as  illustrated  in  figure  2. 

The  results  of  a  parametric  investigation  showed  that  a  performance  giving 
more  than  a  30  percent  increase  in  power  was  achievable.  Furthermore,  fuel 
consumption  could  be  improved  by  10  percent  and  the  CO  emissions  could  be  cut 
to  .1  percent  and  unburnt  hydrocarbons  to  100  ppm,  excellent  by  any  standards 
(Reference  3).  The  precise  reasons  why  this  type  of  configuration  has  such 
results  and  what  the  optimizing  parameters  are  not  understood.  Past  optimiz¬ 
ation  of  the  effect  of  the  balancing  chamber  for  certain  engines  has  been  found 
through  "a  long  and  tedious  process"  of  parametric  variation  (Reference  3). 
Obviously,  this  type  of  guess  work  must  be  eliminated  if  the  concept  of  heat 
balancing  is  to  be  feasible  for  utilization  in  basic  designing  of  engines. 

An  early  step  in  understanding  the  process  associated  with  the  Heat 
Balanced  cycle  was  a  non-steady  flow  analysis  of  the  NAHBE  combustion  process 
completed  in  1978  by  Raghavan  P.  Pandalai.  Using  classical  one-dimensional 
non-steady  flow  solution  techniques,  Pandalai  found  the  exact  time  solution 
of  the  flow  of  a  mass  into  and  out  of  an  idealized  balancing  chamber.  This 
solution  is  shown  in  figure  3.  The  interface  at  non-dimensional  time  zero  is 
simply  the  assumed  separation  point  of  combustible  gases  and  air.  Once 
combustion  is  begun,  a  shock  wave  is  transmitted  into  the  balancing  chamber 
which  accelerates  the  interface  to  motion  towards  the  chamber.  The  shock  wave 
rebounds  from  different  cross-sectional  area  reductions  at  non-dimensional 
distances  equal  to  1.0,  1.78  and  2.2.  The  most  important  point  is  that  the 
interface  is  slowed  and  finally  reversed  in  direction  (region  24  and  31). 

This  indicates  there  might  be  some  type  of  oscillatory  motion  of  the  mass 
volume  into  and  out  of  the  balancing  chamber  (Reference  4). 


It  has  been 


suggested  that  this  oscillatory  motion  could  be  predicted  and  the  exact  non¬ 
steady  solution  for  the  mass  flow  generalized  by  the  Helmholtz  Resonator 
Theorem. 
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THE  HELMHOLTZ  RESONATOR  THEOREM 


A  Helmholtz  Resonator  may  be  visualized  as  a  bottle  with  a  neck  as 
illustrated  in  figure  4.  The  volume  of  fluid  in  the  chamber  is  disturbed  by 
pressure  waves  from  some  source  (imagine  a  small  speaker  humming  at  one 
particular  frequency).  The  fluid  medium  is  alternately  compressed  and  expanded 
by  the  pressure  waves  as  they  move  into  the  opening.  The  fluid  in  the  opening 
may  be  considered  to  act  as  a  mass  element,  a  piston.  There  is  a  stiffness 
related  to  the  size  of  the  opening  and  a  resistance  to  motion  associated  with 
the  radiation  of  energy  in  the  form  of  sound.  Other  resistive  forces,  such 
as  viscosity,  will  not  be  treated  in  this  discussion.  This  system  may  be 
modeled  by  an  analogous  mechanical  system  of  the  type  shown  in  figure  5.  The 
equation  describing  the  displacement  of  the  mass  of  such  a  forced  mechanical 
oscillator  is 

(M)z  +  (B)i  +  (K)z  =  f  (1) 

Now  the  mass,  M,  stiffness,  K,  and  damping,  B,  must  be  described  in  terms  of 
the  acoustic  system. 

The  mass  of  the  mechanical  system  is  transposed  in  acoustic  terms  to 
pQ  SI',  where  pQ  is  the  density  of  the  gas,  S,  the  cross-sectional  area  of  the 
opening,  and  T  the  effective  length  of  the  opening.  This  effective  length 
is  greater  than  the  physical  length  of  the  resonator,  1,  and  related  by  the 
diameter  of  the  opening,  2a,  through  the  equation 

1 1  =1  +  a1  (2) 

where  Al  =  .85a  for  wide  flanged  constrictions  and  Al  =  .6a  for  unflanged 
openings.  The  spring  constant,  K  in  acoustic  terms  is  pQ  c2S2/V  where  c  is 
the  speed  of  sound  through  the  gas  and  V  is  the  volume  of  the  chamber.  The 
damping  constant  or  effective  resistance,  B,  is  found  to  be  pQ  ck2S2/2n  in 
terms  of  acoustic  variables  k,  the  ratio  of  specific  heats  for  the  gas,  and 
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the  others  already  defined.  Finally,  the  driving  force,  f,  is  produced  by 
the  instantaneous  acoustic  pressure  wave  of  amplitude  P  as  it  impinges  on 
the  cross-sectional  area  S.  This  term  is  SP  cos(wt).  Substituting  terms, 
the  differential  equation  describing  the  inward  displacement  of  the  fluid,  y, 

+  <3> 

Just  as  the  acoustic  system  was  modeled  by  a  mechanical  system,  the 
mechanical  system  is  easily  modeled  by  an  electrical  system  such  as  the  one 
in  figure  6.  In  this  case,  the  mass,  M,  or  the  volume  velocity,  u,  for  an 
acoustic  system  which  is  found  through  the  volume  displacement,  X,  by  the 
relationship  u  =  is  analogous  to  the  current  of  the  electrical  system. 

The  voltage  corresponds  to  the  mechanical  and  acoustic  systems  by  V  =  P  cos(wt). 
The  equation  describing  the  electrical  system  made  up  on  an  inductive  element, 
L,  a  resistive  element,  and  a  capacitive  element  of  magnitudes  M,  R,  C 


respectively  is 


TF  *  M  +  h  1dt  *  v 


Again,  the  problem  reduces  to  solving  the  electrical  constants  in  terms  of  the 
acoustic  constants. 

In  order  to  do  this,  inspection  of  the  equation  describing  the  mechanical 
system  is  in  order.  Using  the  equations 

X  =  yS  (5) 

$  =  (6> 


equation  3  reduces  to 


j2v  Pn  Ck2S  Jy  p A  c2S 

1 '  liF  +  —  Ht  +  ~1 —  x  =  SP  C0SM 


7 


and  further  to 


por  d2X  ,  po  ck2  dX  ,  po  °2 
S  dt7  2ir  Ht  V 


X  - 


P  COS(<i)t) 


(8) 


Finally,  since 


and 


1  =  u  =  $T 


(9) 


V  =  P  cos(wt)  (10) 

then  relating  terms  in  equations  4  and  8  it  is  obvious  that 

M  =  p017S  (11) 

R  =  P()  ck2/2Tt  (12) 

J*P0C=/V  (13) 

Thus  the  acoustic  impedance  of  a  Helmholtz  Resonator  may  be  found  through  the 
electrical  analogy  to  be 

Z  =  R  +  j(u)M  -  ^)  (14) 


Equation  9  points  to  the  fact  that  there  is  some  frequency,  w0,  that 
results  in  the  imaginary  term  of  the  acoustic  impedance  to  be  zero  and  resonance 
to  occur,  hence  the  name  Helmholtz  resonator.  This  frequency  is  found  easily 

t>y» 

“o"  -  dr  ■  0  <>5> 


then 


=  (1/MC) 1/2  =  c(S/l • V) 1/2  (16) 

The  similarity  to  the  electrical  system  may  be  carried  even  further  to  the 
actual  frequency  response  of  the  system.  That  is  the  quality  and  sharpness  of 
resonance  may  be  predicted  by  the  analogy  to  the  electrical  Q  point.  That  is 


Q 


o)0M 

nr 


2tt(1'3V/S3)i/2 


(17) 
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It  is  vital  to  note  three  important  points  about  the  preceding  discussion 
(Reference  5). 

Firstly,  the  only  easily  adjustable  variables  of  the  resonator  are  those 
associated  with  its  geometry.  That  is  the  gap  diameter,  the  volume,  cross- 
sectional  area  and  the  length  of  the  entrance  to  the  chamber.  The  theory  is 
completely  general  as  to  the  shape  of  the  volume  except  as  to  a  flanged  or 
unflanged  entrance.  Secondly,  both  the  resonant  frequency  and  quality  of 
response  depend  only  on  geometric  variables.  Finally,  this  theory  generalizes 
the  exact  solution  and  gives  an  average  of  the  fluid  displacement  and  velocity. 
Thus,  such  tedious  methods  of  solution  as  noted  earlier  may  be  side  stepped. 

The  applications  of  the  Helmholtz  Theory  to  the  NAHBE  should  be  apparent 
upon  comparison  of  the  balancing  chamber  in  figure  2  and  Helmholtz  resonator 
shown  in  figure  4.  The  shock  wave  produced  by  combustion  is  a  pressure 
differential  wave  and  can  be  considered  the  driving  wave  as  it  rebound  back 
and  forth  from  the  balancing  chamber  to  the  top  of  the  combustion  chamber  and 
back  into  the  balancing  chamber  (Reference  6). 
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PROBLEM  STATEMENT  AND  METHOD  OF  SOLUTION 


If  the  Helmholtz  Resonator  Theory  is  to  be  accepted  as  a  model  for  the 
NAHBE,  it  is  necessary  to  substantiate  it.  In  order  to  accomplish  the  goal 
of  correlating  the  Helmholtz  Resonator  Theory  and  the  actual  combustion  cycle 
of  the  NAHBE,  it  is  primarily  necessary  to  find  some  way  to  solve  the  Helmholtz 
equation  for  the  volume  displacement  and  velocity.  Once  this  step  is  accomplished, 
the  response  of  the  Helmholtz  system  with  parameters  of  various  balancing 
chamber  geometries  can  be  found.  This  response  can  then  be  scrutinized  to 
see  if  the  Helmholtz  Theory  is  a  useful  and  realistic  model.  The  resonant 
frequency  can  be  readily  compared  to  frequencies  of  shock  waves  moving  in  the 
combustion  chamber. 

Secondly,  physical  data,  either  quantitative  or  qualitative,  of  the  actual 
combustion  process  can  prove  invaluable.  Such  data  could  substantiate  or 
destroy  the  correlation  of  the  Helmholtz  Theory  with  reality.  Especially 
useful  would  be  photographs  of  the  combustion  process.  Such  photographs  might 
capture  the  oscillations  and  the  physical  data  could  be  reduced  to  give  both 
mass  displacement  and  velocities  of  gases  and  combustion  interfaces. 
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COMPUTER  ANALYSIS 


The  solution  of  the  Helmholtz  Equation,  if  it  was  to  be  practical,  must 
be  of  such  a  form  as  to  be  useful  to  some  form  of  computer.  A  usual  method 
of  solving  differential  equations  is  the  LaPlace  transformation  technique. 
Unfortunately,  the  LaPlace  transformation  of  the  Helmholtz  Equation  was  not 
convenient  for  evaluation  on  digital  computers  (Reference  6}.  Another 
solution  to  second  order  differential  equations  that  is  readily  adoptable  to 
digital  computers  is  the  Runge-Kutta  integration  method. 

This  method  was  used  in  the  computer  analysis  of  the  Helmholtz  Resonators 
of  various  design  parameters.  The  digital  computer  analysis  was  accomplished 
by  the  program  "SOLVE"  and  its  subprograms  "EQMOT"  and  "VALUES".  These 
programs,  written  by  Ens.  J.  Brastauskas  and  Midshipman  W.  Johnson,  may  be 
found  in  Appendices  A,  B,  and  C. 

The  Runge-Kutta  method  may  be  used  to  determine  the  solution  of  any  order 
ordinary  differential  equation  of  the  form 


=  f  Hx’  •' 


jn-1 


dx' 


dx1 


n-1 


which  is  comparable  to  the  equation  that  results  from  the  Helmholtz  Theory. 
The  solution  begins  by  transforming  the  variable  so  that  only  first  order 
equations  are  obtained.  This  is  done  by  letting 


Further,  define 


Fi  =  dx 


1  _  d 


p  _  2  _  d£y  _  y 

F2  "  Ux"  ■  HP  -  Y3 

Fn  =  =  f(x,Y,,Y2...Yrl) 

The  Runge-Kutta  method  allows  a  numerical  solution  of  the  differential  equation 
in  the  form  y  (x)  by  approximating  this  solution  over  an  interval  of  x  by 
beginning  at  a  specified  known  starting  value  for  y  (x)  and  using  the  Runge- 
Kutta  recurrence  formula.  Thus,  given  a  second  order  differential  equation, 
such  as  the  Helmholtz  Equation, 


^  Sr> 

Vj  =  y 


Fl  ■  h 


h  ■  Fix,  Y1(  y2)  -0 

d¥, 

which  gives  two  ordinary  coupled  differential  equations  Fj  =  -gj-  and 

F2  -  UP 

For  the  Helmholtz  Equation,  Fj  =  g£ 

F2  =  d2x  *  (-r/m)  Y2  -  (1/cm)  Yx  +  (p/m)  cosu>t. 

The  subprogram  "EQMOT"  gives  values  for  Fj  and  F2  for  the  given  initial 
conditions  and  constants. 
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The  solution  at  particular  points  after  the  initial  conditions,  x  +  ax, 
are  found  by  the  Runge-Kutta  equations  for  a  fourth-order  differential. 

They  are: 

Yli  +  +  2K2  +  2K3  + 

Yi,i+1  § 

\  55  fCx^  Y^,  Y2i)  x 

AX  Ki  \ 

K2  =  f(Xi  +  y,  Y^  +  y,  Y2^  +  y)  AX 

K3  =  f^xi  +  T»  Yii  +  T*  Y2i  +  Ax 


K4  =  f(Xi  +  x,  Yxi  +  K3,  Y2i  +  13)  ax 

Y2i  +1l+2^2+213+14 
Y2,i+1  g 

lj  =  gCx^  Yir  Y2i)  AX 

K  1 

12  =  g(x.j  +  y»  Yn-  +  y,  Yzi  +  y)  AX 

x  K2  1 

13  =  g(Xi  +  £,  y li  +  r  +  T}  x 

14  =  g(x.  +  AX,  Yj.  +  K^,  ■£ '  +  lg)  AX 

where  f(x,  Y1#  Y2j  =  =  Fj  and  g(x,  Yj,  Y2)  =  3^  =  F2 

Thus,  once  an  increment,  H,  is  computed  in  line  420  of  "SOLVE",  the  process 
is  simply  a  loop,  finding  each  K  constant  and  1  constant  (lines  500-620)  and 
then  finding  the  values  of  y  and  at  the  new  point  t  +  H  by  the  formulas 
given  above  (lines  630-640,  Reference  7). 

The  subprogram  "VALUES"  allow^the  user  to  input  geometric  parameters 
of  a  given  engine.  The  values  for  R,  M  and  1/C  are  found  using  thermodynami 
principles  and  plane  geometry  formulas. 


A  plot  of  the  final  solution  for  either  displacement  or  velocity  versus 
time  is  displayed  by  the  library  subprogram  "TEKSUB3".  From  these  plots  the 
maximum  of  the  third  oscillation  is  taken  and  the  value  of  the  frequency  and 
maximum  displacement  was  put  into  data  file  "DATA3"  by  way  of  the  program 
"REVAMP"  written  by  Midshipman  Johnson  (APPENDIX  D).  The  data  file  is  then 
plotted  by  the  library  program  TEKGRAF3  and  curve  fitted  to  give  easily 
readable  data. 


PHYSICAL  ANALYSIS 


In  order  to  obtain  data  of  the  combustion  cycle,  a  time  dependent  film 
analysis  utilizing  Schlieren  photography  techniques  was  undertaken.  In  this 
analysis,  a  Schlieren  apparatus  utilizing  the  Toepler  method  was  used  to 
photograph  density  gradients  occurring  during  combustion  in  a  heat  balanced 
engine.  Such  methods  of  high  speed  photography  in  conjunction  with  Schlieren 
techniques  have  only  been  introduced  very  recently  (Reference  8).  A  short 
review  of  Schlieren  methods  is  in  order. 

A  w*11  known  fact  of  physics  is  that  a  ray  of  light  will  be  deflected 
fro«  'ts  normal  path  as  it  passes  through  a  transparent  medium  of  different 
index  of  refraction  than  that  of  the  medium  in  which  it  was  originally 
traveling.  The  index  of  refraction  is  proportional  to  the  density  of  a  medium, 
and  the  magnitude  and  direction  of  deflection  is  dependent  upon  the  relative 
magnitudes  of  the  densities  of  the  medium  entered  and  exited  by  the  ray  of 
light.  These  principles  are  the  basis  of  the  Schlieren  method.  A  usual 
Schlieren  apparatus  is  illustrated  in  figure  7. 

A  point  light  source  at  S,  is  stationed  at  a  distance  equal  to  the 
focal  distance,  f,  of  a  concave  mirror,  Mx.  This  results  in  a  line  of 
parallel  ray  reflected  off  the  mirror  and  through  the  working  section.  Any 
disturbances  to  the  density  of  the  medium  in  the  section  will  cause  a 
deflection  of  the  rays  (dotted  lines).  Any  rays  that  are  undisturbed  are 
passed  through  in  a  parallel  and  straight  course.  Both  sets  of  rays  are 
reflected  by  another  concave  mirror,  M2.  The  two  sets  of  rays  fall  on  the 
screen  Q.  The  Toepler  method  utilizes  a  knife  edge  at  the  focal  point  of 
the  second  mirror.  The  knife  edge  effectively  blocks  out  any  rays  that  are 
bent  too  far  thus  creating  a  shadow.  Any  light  rays  bent  in  the  opposite 
direction  cause  an  addition  of  light  on  the  screen  resulting  in  a  light  spot. 
Thus,  any  density  gradients  of  the  medium  may  be  viewed  (Reference  9). 
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Such  an  apparatus  was  constructed  for  photographing  the  combustion  of 
the  NAHBE.  A  square  engine  with  two  glass  walls  (fused  silicon)  enabled  the 
light  rays  to  pass  through  the  section  and  not  be  effected  by  curved  glass 
walls.  The  engine  was  run  at  1,000  rpm  in  spark  ignition  mode  and  was  the 
test  section.  The  light  source  was  a  500  watt  xenon  arc  lamp  focused  through 
a  converging  lens.  A  small  slit,  acting  as  a  horizontal  polarizer,  was  placed 
just  behind  the  focal  point  of  this  lens  so  as  to  eliminate  any  chromatic 
aberation.  The  actual  set  up  is  illustrated  in  figure  8  and  the  test  section 
(square  engine)  is  shown  in  figure  9. 

In  order  to  capture  the  combustion  sequence  and  slow  it  down  so  that 
useful  data  might  be  observed,  a  high  speed,  Fastek  motion  picture  camera  was 
used.  Previous  studies  showed  that  a  camera  speed  of  2,000  frames  per  second, 
and  a  corresponding  shutter  speed  of  1/6, 000th  of  a  second  were  best.  This 
allows  a  picture  to  be  taken  about  every  .5  milliseconds.  A  light  meter 
reading  of  about  15.5  is  optimal  also  (Reference  8). 


RESULTS 


Computer  Analysis 

The  generalized  time  response  of  several  Helmholtz  resonators  having 
geometries  corresponding  to  several  NAHBE  caps  now  in  use  follows.  Figures 
10-21  show  the  time  response  for  frequencies  ranging  from  3000  to  40,000 
radians  per  second  of  the  square  engine  with  a  NAHBE  cap  having  dimensions 
stated  in  table  II. 


Table  II 

NAHBE  Square  Engine  Parameters 


Compression  ratio 
cap  thickness 
gap 

cross-sectional  area 
volume 

resonant  frequency 
quality  number 


y  =  3.5 
1  =  .1  inch 
2a  =  .0625  inch 
S  =  .101  sq.  in. 
V  *  .393  cu.in. 
u0  =  23919.2 
Q  =  6.99 


The  pressure  of  the  wave  was  the  only  problem  in  the  analysis.  The 
value  for  the  pressure  of  an  acoustic  wave  can  be  stated  in  various  ways. 
Therefore,  a  unit  magnitude  value  of  1  acoustic  pressure  was  used  for  the 
pressure  value  necessary  for  the  solution.  This  only  changes  the  magnitude 
of  the  driver  function  which,  in  return,  effects  the  magnitude  of  the  displace¬ 
ment.  This  is  of  no  consequence,  since  it  is  the  shape  of  the  time  response 
solution  desired  and  not  the  magnitude.  In  any  case,  these  magnitudes  may 
be  considered  as  relative. 

After  performing  similar  sets  of  solutions  for  various  caps,  the  maximum 
displacement  occurring  during  the  third  cycle  was  plotted  against  the  frequency. 
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These  plots  are  found  in  figures  22  through  25.  It  should  be  noted  that  the 
plot  fits  were  used  only  to  give  a  general  idea  of  the  shape  of  the  curve. 

A  close  inspection  of  these  plots  show  interesting  results.  The  most  obvious 
difference  in  the  plots  is  the  frequency  at  which  the  peak  displacement  occurs. 
Further  scrutinization  shows  that  the  peak  becomes  sharper  and  slopes  more 
steeply  for  each  cap.  Lastly,  it  is  evident  that  the  displacement  before  the 
peak  is  significantly  greater  than  those  after.  This  trend  is  especially  true 
for  the  smaller,  squatter  curves. 

Table  III  shows  the  computed  values  of  resonant  frequency,  quality  number 
for  balancing  chambers  whose  cap  thicknesses  were  varied.  The  volume  and 
cross-sectional  area  remained  constant  for  all  the  runs. 


Table  III 


1 

a 

% 

.1" 

6.99 

23919.2  rad/sec 

.2" 

14.2 

18199.7  rad/sec 

.3" 

24.1 

15267.2  rad/sec 

It  is  evident,  in  consideration  of  the  data  in  table  III,  that  the  peaks  of 
the  frequency  plots  occur  at  resonance  for  the  chambers.  This  explains  the 
shift  in  the  curves  as  being  due  to  a  change  in  the  resonant  frequency  due  to 
the  variation  or  the  cap  thickness.  Further  investigation  shows  that  the 
quality  number,  Q,  gives  an  idea  as  to  the  intensity  of  the  response  at 
resonance.  A  relatively  low  Q  for  the  1  =  .1"  geometry  shows  that  the 
frequency  width  of  the  response  at  .707  times  the  maximum  (around  .0000277) 
is  approximately  from  18000  rad/sec  to  27000  rad/sec  or  a  loosely  termed 
bandwidth  of  9000  rad/sec.  In  contrast,  the  cap  of  thickness  .3"  and  Q  *  24 
has  a  bandwidth  of  5000  rad/ sec. 
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The  Cooperative  Fuel  Research  (CFR)  engine  has  been  used  in  many  para¬ 
metric  variation  experiments.  The  cap  used  in  the  CFR  has  a  geometry  of 
1  =  .0749",  S  =  .7967  sq.  in.  and  V  =  2.398  cu.  in.  resulting  in  =  33937.8 
and  Q  =  .5895.  The  same  relationships  discussed  above  are  concurred  by  the 
CFR  data.  The  bandwidth  is  found  to  be  14500  rad/sec  approximately. 

Physical  Analysis 

The  photographic  results  of  the  heat  balanced  cycle  in  the  square  engine 
proved  hard  to  reduce  to  quantitative  data.  The  subjective  points  of  the  film 
are  easily  seen  when  viewed  in  motion.  Unfortunately,  limited  space  allows 
only  one  combustion  sequence  of  one  of  the  films  taken  to  be  discussed.  Eight 
other  films  of  excellent  quality  were  taken.  Such  parameters  as  volume  and 
flanging  were  varied.  The  results  discussed  are  typical  of  all  the  film 
analyses  done  on  the  square  engine. 

It  should  be  noted  that  the  photographs  appearing  in  figures  26  through 
64  are  the  negatives  of  color  prints.  Therefore,  the  light  portions  are  dark 
and  the  dark  portions  light.  Furthermore,  it  must  be  kept  in  mind  that  these 
pictures  are  blown-up  from  16mm  film.  As  a  result,  the  quality  of  the  picture 
has  been  reduced.  Fortunately,  enough  detail  exists  for  a  reasonable  discussion. 

The  cap  used  in  these  films  was  square  in  cross-section.  The  thickness 
of  one  flange  was  approximately  .1  in.  and  the  other  side  was  .2  in.  The 
volume  and  cross-sectional  area  of  the  gap  were  kept  constant  at  .393  cu.  in. 
and  .101  sq.  in.  respectively.  This  configuration  allowed  testing  and  filming 
of  two  different  geometries  under  identical  conditions.  Furthermore,  the 
geometries  are  the  same  as  those  of  the  computer  analysis  models. 

The  sequence  begins  as  spark  ignition  takes  place  in  the  chamber.  The 
flame  front  can  be  seen  to  expand  through  figure  36.  In  figure  35  and  36, 
shadows  resembling  ripples  in  a  pond  can  be  seen  under  the  left  side  of  the 


cap.  These  rings  are  most  likely  the  reflected  shock  wave.  Unfortunately, 


r 


the  speea  of  the  shock  wave  is  so  fast  that  it  is  impossible  to  tell  whether 
this  is  the  first,  second  or  third  reflection.  The  initial  combustion 
continues  to  its  completion  until  about  figure  50  where  it  looks  as  though 
the  combustion  is  completed.  In  figure  51,  the  rings  are  again  seen,  but  now 
above  the  cap  on  the  left  side.  This  reaction,  which  most  likely  is  the 
oscillating  of  the  mass  of  air  in  the  gap,  increases  until  figure  55.  The 
combustion  shows  regeneration  and  spreads  across  the  cap  in  figures  53  through 
58.  At  approximately  figure  58  the  secondary  combustion  of  the  left  side 
meets  an  impass.  This  is  a  result  of  a  very  strong  reaction  occurring  from 
the  right  chamber  at  this  time.  In  figure  62  the  secondary  combustion  is  all 
but  completed  and  is  finally  quenched  as  the  exhaust  valve  begins  to  open  in 
figures  63  through  64. 

From  a  subjective  point  of  view,  the  film  does  not  necessarily  substantiate 
the  use  of  the  Helmholtz  model.  It  does  show  two  distinct  combustion  processes. 
Also,  it  shows  there  is  some  sort  of  flow  from  under  the  cap  causing  the 
secondary  reaction.  A  closer  evaluation,  however,  shows  that  the  film 
analysis  can  be  used  for  quantitative  data  of  some  kind. 

It  is  known  that  the  combustion  cylinder  is  charged  with  a  propane-air 
gas  solution  at  approximately  bottom  dead  center.  The  volume  of  the  cylinder 
at  this  point  is  found  by  the  product  of  the  cross-sectional  area  and  the 
length  of  the  cylinder  (from  top  of  the  cap  to  top  of  the  cylinder).  It  is 
also  known  that  the  compression  ratio  (the  ratio  of  volume  at  bottom  dead 
center  and  top  dead  center)  is  3.5.  Thus,  the  volume  of  gas  when  the  piston 
is  at  top  dead  center  is  the  product  of  the  cross-sectional  area  and  the 
length  of  the  combustion  chamber,  H.  Mathematically  stated,  this  is 

VBDC  =  HBDC  xA 

VTDC  =  HTDC  xA 
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The  compression  ratio  then  is 


VBDC  _  ,  -  _  HBDC 

W  ~  u 

VTDC  "TDC 


(20) 


The  volume  ratio  is  also  associated  with  the  temperature  of  the  gases  which 
undergo  isentropic  compression  or  expansion  by  the  thermodynamic  formula 
(Reference  10). 

T2  =  Tx  CV1/V2)k"1  (21) 

All  temperatures  are  in  absolute  degrees  and  k  is  again  the  ratio  of  specific 
heats  of  the  gas  undergoing  compression  or  expansion.  For  the  propane-air 
mixture,  it  is  considered  equal  to  the  value  for  air  of  1.4  (Reference  11). 

Therefore,  given  an  initial  temperature  at  the  beginning  of  compression, 
the  temperature  of  the  gas  may  be  found  at  any  other  time  if  the  volume  or 
length  ratio  is  known.  If  this  procedure  is  followed  for  the  square  engine, 
and  assuming  an  initial  temperature  of  76°  or  536°R,  the  temperature  at  the 
end  of  compression  should  be  about  884°R.  Actually,  losses  due  to  non- 
reversible  processes  will  cause  it  to  be  a  little  higher.  It  is  at  this  point 
that  combustion  begins. 

Previous  studies  have  shown  a  wave  of  Mach  1.26  is  generated  by  the 
combustion  gases  (Reference  4).  To  be  conservative,  let  us  assume  the  average 
value  of  all  waves  to  be  a  Mach  of  1.1.  The  speed  of  sound  in  any  gas  is 
associated  with  the  temperature  of  the  gas  by  the  formula  (Reference  12). 

a  =  (kRgT)1/2  (22) 

where  R  is  the  gas  constant  for  any  gas  and  g  the  acceleration  constant  of 
gravity.  For  Old  English  units  and  air  the  equation  reduces  to 

a  =  49.01  (T)1/2  (23) 

Knowing  the  temperature  of  the  gas  then,  the  speed  of  the  wave  in  feet  per 
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second  may  be  determined.  It  is  then  simple  ari thematic  to  find  the  time 
required  for  the  wave  to  begin  from  the  top  of  the  cylinder,  reflect  from  the 
cap  and  return  to  the  cylinder  top,  knowing  the  distance  to  the  cap.  That  is 
the  time,  t,  required  to  complete  this  one  oscillation  is 

t  -  1-1  02  in/ft) (a  ft/sec) 
t  Hlru 

The  frequency,  in  radians  per  second,  of  this  oscillation  is  found  by 

_  2 it  rad/cycle 

“  '  ~  sec /cycle  (25) 

Thus,  the  position  of  the  piston  can  be  used  to  find  the  frequency  of 
the  oscillating  wave  occurring  in  the  cylinder.  This  calculation  was  done, 
for  the  square  engine,  by  the  computer  program  "FREQ"  (APPENDIX  E).  The 
resulting  data  was  placed  into  a  data  file  "FREQDATA"  and  plotted  (figures  65 
and  66). 

When  these  frequencies  are  compared  to  those  found  by  the  computer 
analysis, (figures  22  through  24)  they  are  low  and  should  not  cause  any 
significant  reaction  except  possibly  at  the  beginning  of  the  combustion. 

Closer  consideration  of  the  problem,  however,  shows  a  shortcoming  in  the 
previous  analysis.  The  temperature  at  the  end  of  compression  was  considered 
to  stay  constant  throughout  the  rest  of  the  cycle,  but  the  volume  of  the 
cylinder  is  continually  increasing,  reducing  the  temperature  by  the  same 
formula  as  before  (equation  21).  This  would  result  in  slower  velocities  and 
even  lower  frequencies.  This  is  evident  by  the  results  of  the  revised  program 
"FREQ"  (APPENDIX  F)  in  figures  67  and  68. 

One  last  failacy  in  the  model  is  apparent.  As  propane  bums  in  the  air, 
the  temperature  of  the  gas  in  the  cylinder  will  increase.  The  temperature  of 
combustion  for  propane  is  between  1374°R  and  1518°R.  An  average  value  of 
1446°R  is  assumed.  It  is  obvious  that  the  metal  of  the  cylinder  and  the 
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cooling  system  of  the  engine  will  reduce  this  temperature  significantly.  To 
be  conservative,  again,  it  is  assumed  that  the  temperature  is  increased  by 
20  percent  of  the  temperature  of  combustion,  that  is  the  temperature  in  the 
cylinder  is  raised  by  290°R  throughout  the  process. 

With  this  last  revision  of  "FREQ"  (APPENDIX  6),  the  output  is  shown  in 
figure  69.  When  the  values  of  resonant  frequencies  for  the  two  geometries 
studied  in  the  film  analysis  and  the  chamber  lengths  for  these  frequencies  are 
compared,  interesting  results  are  found.  Such  a  comparison  is  compiled  in 
table  IV. 

Table  IV 

1  ^o  H 

.1"  23919.2  rad/sec  1.5" 

.2"  18199.7  rad/sec  1.9" 

Using  bottom  dead  center  as  a  guiding  distance,  the  photographs  may  be 
scaled  no  matter  how  large  or  small  they  are  printed.  Thus,  a  chamber  length 
of  1.5  inches  is  found  in  figure  54  of  the  film  sequence.  The  chamber  length 
of  1.9  inches  is  found  at  around  figure  58.  These  finds  of  the  computer  model 
in  "FREQ"  and  the  subjective  data  of  the  film  analysis  are  corroborative. 
Similar  studies  on  past  film  analysis  obtain  comparable  results  (Reference  8). 


SUMMARY  OF  RESULTS 


Computer  analysis  shows  that  the  square  engine  caps  modeled  by  the 
Helmholtz  Theory  had  a  correlation  between  geometry  and  resonant  frequency, 
quality  of  reaction,  and  bandwidth  or  reaction. 

Physical  analysis  by  high  speed  filming  shows  a  definite  correlation 
between  the  computer  analysis  prediction  and  actual  reaction.  The  film 
captured  every  aspect  of  the  combustion  process  and  parts  of  the  wave  inter¬ 
action. 
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CONCLUSION 


It  is  evident,  by  the  correlation  of  the  computer  and  physical  data, 
that  the  Helmholtz  Resonator  Theory  describes  the  general  time  dependent 
reaction  of  the  NAHBE  cap.  It  may  be  possible  to  predict  the  actual  physical 
reaction  occurring  during  combustion  using  the  Helmholtz  Theory.  Furthermore, 
the  Helmholtz  Theory  may  be  used  to  optimize  the  reaction  of  heat  balanced 
processes  now  in  use  in  the  NAHBE.  It  may  also  possibly  be  used  as  a  design 
tool  to  create  optimum  output  given  operating  specifications  of  an  engine. 

Further  studies  should  be  made  to  find  the  relationship  between  the 
intensity  of  the  reaction  and  output  power.  That  is  to  say,  the  study  should 
find  which  is  more  desirable,  a  longer  weaker  reaction  corresponding  to  a 
low  Q  value,  or  a  strong,  short  reaction  typical  of  high  Q  values. 

The  film  analysis  of  combustion  gives  valuable  insight  into  the  actual 
process.  Quantitative,  as  well  as  qualitative,  data  may  be  obtained  and  used 
in  future  studies.  Reasons  for  such  side  effects  as  low  hydrocarbon  emissions 
are  evident  by  viewing  the  extremely  long  combustion  of  the  gases.  Film 
analysis  has  been  refined  to  the  point  that  it  can  easily  be  used  for  any  type 
of  future  research. 
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1O0  REM  THIS  PROGRAM  WILL  ALLOW  YOU  TO  STORE  THE  W  VALUES  AMD  THE 
110  REM  CORRESPONDING  MAXIMUM  DISPLACEMENT  (3RD  CYCLE)  INTO  A 
120  REM  DATA  FILE,  "DAT A3".  THIS  THEN  COULD  EE  USED  WITH 
130  REM  OLD  L.  IG***: TEKGRAF3  TO  PLOT  THE  FREQUENCY  RESPONSE  FOR 


140  REM  THE  CERTAIN  GEOMETRY 
130  TILE  #1:  "DATA3" 

ISO  PRINT  "INPUT  W" 

170  INPUT  W 


130  PRINT  "INPUT  MAXIMUM  DISPLACEMENT 
ISO  INPUT  D 
00  PRINT  #1:  W;  ",  " ;  D 


( 3RD  CYCLE ) " 


10  END 
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RAGE 


100 

110 

120 
1 30 
140 
150 
160 
170 

1  SO 
1 90 
200 
210 
220 
230 
240 
250 
260 
270 

2  SO 
290 
300 
310 
320 

330' 

340 
350 
'  360 
370 


REM  THIS  PROGRAM  IS  DESIGNED  TO  FIND  THE  FREQUENCY  THAT  A  MACH  WAVE 
REM  WILL  OBTAIN  IF  PUT  IN  A  COMBUSTION  CHAMBER  OF  A  CERTAIN  VOLUME 
REM  COMPARED  TO  THE  VOLUME  AT  THE  TIME  OF  COMBUSTION. 

REM  IT  WILL  BE  NECESSARY  TO  INPUT  THE  LENGTH  OF  THE  COMBUSTION 
REM  CHAMBER  AT  THE  TIME  THE  FREQUENCY  IS  SOUGHT  FOR.  THE  INITIAL 
REM  CONDITIONS  FOR  THE  CHAMBER  AT  THE  BEGINNING  OF  COMPRESSION  ARE 
REM  T0=536R<  P0=14.  7PSI.  THE  APPROX  I MAT I ON  THAT  THE  WAVE  IS 
REM  TRAVELING  AT  1.  1  TIMES  THE  SPEED  OF  SOUND  THAT  OCCURS  IN  THE 
REM  CHAMBER  IS  EVIDENT.  OTHER  ASSUMPTIONS  ARE:  ALL  GASES  I-N  THE 
REM  COMBUSTION  CHAMBER  ARE  OF  THE  SAME  CHARACTERISTICS  AS  AIR. 

F I LE# 1 : " FREQDAT A "  *  * 

SCRATCH# 1 

PRINT" INPUT  ALL  DISTANCES  IN  INCHES” 

PRINT" INPUT  THE  COMPRESSION  RATIO" 

INPUT  V 

PRINT  "INPUT  THE  BOTTOM  DEAD  CENTER  CHAMBER  LENGTH" 

INPUT  C 

LET  T2=(536*V)  '\  4 

PRINT  "INPUT  THE  TOP  DEAD  CENTER  CHAMBER  LENGTH" 

INPUT  LI 

FOR  L=L 1  TO  C  STEP  .  1 
LET  A-49.  01*12*  <  <  T2 )  "\  5 ) 

LET  !J=1.  1*A 

LET  R=L/U 

LET  W-3.  1415/R 

PRINT  "W  IS"; W;  " RAD/ SEC" 

FEIMT#1 :  Wi  "»  '*;  L 
NEXT  L 
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REM  THIS  PROGRAM  IS  DESIGNED  TO  FIND  THE  FREQUENCY  THAT  A  MACH  WAVE 
REM  WILL  OBTAIN  IF  PUT  IN  A  COMBUSTION  CHAMBER  OF  A  CERTAIN  VOLUME 
REM  COMPARED  TO  THE  VOLUME  AT  THE  TIME  OF  COMBUSTION. 

REM  IT  WILL  BE  NECESSARY  TO  INPUT  THE  LENGTH  OF  THE  COMBUSTION 
REM  CHAMBER  AT  THE  TIME  THE  FREQUENCY  IS  SOUGHT  FOR.  THE  INITIAL 
REM  CONDITIONS  FOR  THE  CHAMBER  AT  THE  BEGINNING  OF  COMPRESSION  ARE 
REM  T0=5S6Ri  P0=14.  TPS I.  THE  APPROXIMATION  THAT  THE  WAVE  IS 
REM  TRAVELING  AT  1.  1  TIMES  THE  SPEED  OF  SOUND  THAT.  OCCURS  IN  THE 
REM  CHAMBER  IS  EVIDENT.  OTHER  ASSUMPTIONS  ARE:  aLl  GASES  IN  THE 
REM  COMBUSTION  CHAMBER  ARE  OF  THE  SAME  CHARACTERISTICS  AS  A'lR. 
FILEttl: "FREQDATA"  .  - 

SCRATCH# 1  * 

PRINT" INPUT  ALL  DISTANCES  IN  INCHES" 

PRINT" INPUT  THE  COMPRESSION  RATIO" 

INPUT  V 

PRINT  "INPUT  THE  BOTTOM  DEAD  CENTER  CHAMBER  LENGTH" 

INPUT  C 

LET  T2=  <  S36*V )  "\  4 

PRINT  "INPUT  THE- TOP  DEAD  CENTER  CHAMBER  LENGTH" 

INPUT  LI 

FOR  L=L1  TO  C  STEP  .  1 
LET  TS=<  T2XL1/L)  "/  4) 

LET  A=4S.  01*12*!  !T3)'\  5) 

LET  U-l.  1*A 

LET  R-L/’J 

LET  W=3.  141S/R 

PRINT  "W  IS"; W, "RAD/SEC" 

PRINTttl:  W;  ">  "j  L 

NEXT  L 

END 


Appendix 


xxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxx 


•: i  2 

2  2 

a  2 

22222 
g"  g 


2222 
2  « 

2  e 

2  « 

2  fi 

2222 


2222C 

2 

2222 

2 

2 

2 


xxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxx; 


1.1 


IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

;xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


(3  *2*5  <3  <5 

<2  <2  <2 

aaa 

aaa 

aaa 

a 

a 

ea 

a 

a 

a 

a 

a 

<§ 

a 

<2 

a 

a 

aa 

a  aa 

a 

a 

a 

a  -  a 

a 

ee  1 

a 

a 

a 

a 

a 

a 

a  a 

a  i 

a 

see 

a 

aaa 

aaa 

aaa 

aeea 

aaaaa 

aaa 

s 

a 

a 

a  a 

a 

a 

eaaa 

a  a 

aeaa 

a 

a  a 

\ 

a 

a 

a 

a  aa 

a  a 

aaaaa 

•»* 

\ 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
xxxxxxxxxxxxxxxxxx  xxxxxxxxxxxxxxxxxxxxxxxxxx 


35 


04/01/31 


04/01/31 


10:  42:  51 


PAGE  1 


EG 


•0  REM  THIS  PROGRAM  IS  DESIGNED  TO  FIND  THE  FREQUENCY  THAT  A  MACH  WAVE 

.0  REM  WILL  OBTAIN  IF  PUT  IN  A  COMBUSTION  CHAMBER  OF  A  CERTAIN  VOLUME 

■0  REM  COMPARED  TO  THE  VOLUME  AT  THE  TIME  OF  COMBUSTION.  4 

50  REM  IT  WILL  BE  NECESSARY  TO  INPUT  THE  LENGTH  OF  THE  COMBUSTION 

10  REN  CHAMBER  AT  THE  TIME  THE  FREQUENCY  IS  SOUGHT  FOR.  THE  INITIAL 

50  REM  CONDITIONS  FOR  THE  CHAMBER  AT  THE  BEGINNING  OF  COMPRESSION  ARE  € 

»0  REM  T0=536R»  P0=14.  7PSI.  THE  APPROXIMATION  THAT  THE  WAVE  IS 

70  REM  TRAVELING  AT  1.  1  TIMES  THE  SPEED  OF  SOUND  THAT  OCCURS  IN  THE 

30  REM  CHAMBER  IS  EVIDENT.  OTHER  ASSUMPTIONS  ARE:  ALL  GASES  IN  THE  € 

■■0  REM  COMBUSTION  CHAMBER  ARE  OF  THE  SAME  CHARACTERISTICS  AS  AIR. 

30  FILE#1: "FREQDATA"  *  - 

10  SCRATCH# 1  1 

20  PR I NT "INPUT  ALL  DISTANCES  IN  INCHES" 

30  PRINT" INPUT  THE  COMPRESSION  RATIO" 

40  INPUT  V  C 

50  PRINT  "INPUT  THE  BOTTOM  DEAD  CENTER  CHAMBER  LENGTH" 

SO  INPUT  C 

70  LET  T2=*<536*V)-\  4  f 

SO  PRINT  "INPUT  THE  TOP  DEAD  CENTER  CHAMBER  LENGTH" 

90  INPUT  LI 

00  FOR  L=L1  TO  C  STEP  .  1  f 

10  LET  T3=<T2*<L1/L)'\  4) 

20  LCT  T4=T3+2S9.  2 

30  LET  A=49.  01*12*  <  <T4>-\  5)  '  6 

40  LET  U=l.  1*A 
50  LET  R=L/U 

•SO  LET  W«3.  1415/R  Q 

70  PRINT  "W  IS"; W; "RAD/SEC" 

SO  PR  I  NT#  1:  W;  ",  L 

90  NEXT  L  ,  Q 

00  END 
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